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TEMA 1: Reactores catalíticos HETEROGÉNEOS 


O. Resumen 


El siguiente tema consta de dos partes. La primera es un estudio de las 
reacciones heterogéneas catalizadas por sólidos, donde se incluye un seguimiento 
detallado de las ecuaciones cinéticas de las mismas según cual sea la etapa 
determinante de la velocidad de reacción: difusión externa desde el entorno hasta la 
superficie, adsorción, reacción, desorción, difusión interna a través de una pastilla 
porosa, efectos caloríficos o combinaciones de varias de las anteriores. La segunda 
trata sobre el diseño de los reactores catalíticos, de los cuales existe una gran 
variedad, aunque de todos ellos solo se va a estudiar los reactores de lecho fijo por 
etapas. El estudio se centrará también en los reactores adiabáticos con 


enfriamiento entre etapas. 


1. Introducción 


La velocidad de una reacción puede ser modificada por la presencia de unas 
sustancias, que normalmente no son ni reactivos iniciales ni productos. Estas 
sustancias reciben el nombre de catalizadores, y producen un aumento de la 
velocidad de reacción. Los catalizadores pueden ser sólidos o líquidos. Los 


catalizadores sólidos pueden perder su actividad catalítica con el tiempo. 


Un catalizador puede hacer variar la velocidad de reacción de un proceso en 
miles o millones de veces, por lo que suele ser muy interesante su utilización en 
procesos industriales (producción de ácido sulfúrico, amoníaco, etc). Las 


características de un catalizador se podrían resumir en: 


1. Selectividad de los catalizadores. Se refiere a la capacidad que tienen estas 


sustancias de actuar en ciertas reacciones y de no hacerlo en otras distintas. 


2. Un catalizador aumenta la velocidad de reacción o favorece una reacción 
frente a otra, pero en ningún caso determina el equilibrio o el punto final de la 


misma, que vendrán condicionados siempre por cuestiones termodinámicas. 


3. Para utilizar un catalizador en una reacción determinada, se hace 
necesario llevar a cabo un número elevado de ensayos y pruebas, hasta dar con el 


más adecuado, dentro de una lógica. 
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4, La actividad de un catalizador no sólo depende de su constitución química, 
sino también de su estructura física o cristalina, ya que el catalizador puede perder 
su actividad por encima de ciertas temperaturas. Por lo tanto, un estudio 
exhaustivo de un catalizador debe incluir una investigación de la superficie del 
mismo. Por ello se ha de procurar que los catalizadores tengan una gran superficie 


sólida activa por unidad de volumen. 


5. En presencia de un catalizador las moléculas reaccionantes sufren un 
debilitamiento de sus enlaces o forman productos intermedios en las proximidades 
de la superficie del sólido. Se han propuesto varios modelos que intentan explicar 


este comportamiento: 


a) el producto intermedio es una asociación de la molécula reactivo con una 


región de la superficie del sólido 


b) las moléculas de reactivo se encuentran bajo la influencia de las fuerzas de 


la superficie, moviéndose en un entorno cercano a ella 


c) se forman radicales libres en la superficie del catalizador, que se mueven 
hacia el gas provocando una serie de reacciones, hasta que finalmente se 
destruyen; en esta última teoría, la reacción va a tener lugar en el seno del 


gas, siendo la superficie del catalizador un generador de radicales libres. 


6. Según la teoría del estado de transición, el catalizador reduce la barrera de 
energía potencial que se debe sobrepasar para que los reaccionantes pasen a 
productos. Esto provocará a su vez una disminución de la energía de activación, y 
por lo tanto, un aumento de la velocidad de la misma, como se puede observar en la 


Figura 1. 


Para poder llevar a cabo un diseño de reactores para este tipo de reacciones, 
será necesario estudiar detalladamente la ecuación de velocidad en estos procesos a 


partir del estudio de la cinética de la reacción. 
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Sin catalizador, el 
complejo tiene más energía 


Energía 


Estado inicial 


La barrera de energía 
desciende con catalizador, 


Productos 
permitiendo velocidades mayores P 


Reactivos Complejo Productos 


Transcurso de la reacción 


Figura 1. Representación de la acción de un catalizador. 


2. Ecuación de velocidad 


En la Figura 2, se muestra un esquema de los tipos de reactores ideales. Un 
balance en unidades de flujo (kmol/s) del componente A, y considerando "flujo de 
salida - flujo de entrada + flujo de acumulación = flujo de generación", conduce a 


las siguientes expresiones: 


— — 

==] 

Reactor continuo 
Reactor discontinuo de flujo pistón Reactor continuo 
de tanque agitado de tanque agitado 


Figura 2. Tipos de reactores ideales. 


a) Para el RDTA, aplicado a todo el volumen 
dNA/dt=8gA (2.1) 
siendo Na los kmol de A que hay dentro del reactor, y ga los kmoles generados/s en 
este reactor, y por tanto es una propiedad extensiva. 


b) Para un elemento de volumen del reactor de flujo de pistón: 
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na+dna-na=d8a (2.2) 


siendo na el flujo molar de entrada (kmol/s), na+dna el flujo molar de salida y dga la 


generación de A. 
c) Para el reactor continuo de tanque agitado: 
Ma -Mao = Ba (2.3) 


En sistemas catalíticos, se puede encontrar la velocidad de reacción 


expresada de muchas formas. 

1. Basada en la superficie externa del catalizador 

ra= ga/S [kmoles reaccionados/m2 superficie externa catalizador:s] 
siendo S la superficie externa del catalizador. 

2. Basada en el peso de las partículas del catalizador 

ra=ga/M [kmoles reaccionados/kg de catalizador:s] 
siendo M la masa de catalizador. 

3. Basada en el volumen de la mezcla reactivo 

r',= ga/V [kmoles reaccionados/m* mezcla react.s] 
siendo V el volumen de la mezcla reactivo. 

4. Basada en el volumen de las particulas de catalizador 

r"»= ga/Vo [kmoles reaccionados/m:* part.s] 
siendo V, el volumen de las partículas. 


A partir de las expresiones anteriores, despejando ga, se puede demostrar la 


siguiente equivalencia: 
8,=115=r,W =r,V =r, V, (2.4) 


En catalizadores porosos, la reacción transcurre en toda la interfase sólido- 
fluido tanto en la superficie externa, como dentro de los poros de la partícula. Para 
deducir las expresiones de las velocidades de reacción, es necesario tener en cuenta 
las distintas etapas (tanto fisicas como químicas) que tienen lugar en el proceso de 
reacción y que pueden ofrecer resistencia a la misma, ya que según predomine una 
u otra de estas etapas se tendrán distintas expresiones de la velocidad. En la Figura 
3 se muestra un esquema de las etapas que debe sufrir un reactivo A para llegar a 


la superficie. 
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Figura 3. Etapas de reacción 


Las etapas más lentas que las demás, que limitan el valor de la velocidad de 


reacción se las conoce como controlantes o limitantes. 
Las distintas etapas son las siguientes: 


1. Difusión externa (o resistencia de la película): difusión de los reactivos 


hacia la superficie exterior del catalizador. 


2. Difusión interna (o resistencia a la difusión en los poros): difusión de los 
reactivos hacia el interior del catalizador, donde se va a producir la reacción 


mayoritariamente. 


3. Adsorción química (o física en ciertas ocasiones) en la superficie: durante la 


unión de las moléculas del reactivo en algún lugar de la superficie del catalizador. 


4. Reacción química superficial: reacción en la superficie de la partícula de 


catalizador. 


5. Desorción química en la superficie: durante la separación de las moléculas 


de producto de la superficie del catalizador. 


6. Difusión interna (o resistencia a la difusión en los poros): difusión de los 


productos hacia el exterior del catalizador. 


7. Difusión externa (o resistencia de la película): durante la salida de los 


productos hacia el seno de la fase gaseosa. 


Otros factores que pueden influir 
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8. Flujo de calor en la propia partícula: aparece fundamentalmente en 
reacciones rápidas exotérmicas o endotérmicas, y provoca gradientes de 


temperatura dentro de la partícula. 


9. Flujo externo de calor: provoca un gradiente de temperaturas entre la 
superficie exterior de la partícula y la corriente de fluido exterior. En este caso, el 
catalizador podría tener temperatura uniforme en todos sus puntos, pero distinta a 


la temperatura del seno del fluido exterior. 


De todas las etapas de la reacción, la de difusión interna es algo especial, 
puesto que no puede ser etapa controlante, en el sentido de que si es muy lenta la 
velocidad vendrá dada por la de difusión interna, puesto que la reacción podría 
producirse en la superficie del sólido. Por ello se dice que la resistencia a la difusión 


£- 


interna está "en paralelo" con las otras resistencias, todas ellas "en serie". 


3. Ecuación de velocidad para las etapas químicas como etapas 


limitantes 


3.1. Mecanismos de catálisis 


La naturaleza exacta de las razones por las que se forma fácilmente el 
complejo intermedio, adsorbido sobre la superficie del catalizador, es un tema de 
investigación. Se pueden prever requerimientos estructurales para la formación y 
ruptura de enlaces, y ello ha conducido a la formulación de una serie importante de 
teorías. Como se ha comprobado, la transferencia de electrones está implicada en la 
formación de los enlaces del compuesto intermedio, ya por ello el carácter 


electrónico de los catalizadores es indudablemente importante. 


No se va a pretender dar una explicación satisfactoria, ya que no la existe, de 
los diferentes mecanismos, sino que se van a proponer algunos ejemplos con las 


correspondientes explicaciones que presentan los investigadores en la bibliografia. 


Por ejemplo, los catalizadores ácidos, tales como sílice y alúmina, pueden 
aparentemente actuar como ácidos de Lewis (aceptor de electrones) o ácidos de 
Brónsted (dador de protones). El mecanismo de reacción para el craqueo de un 


hidrocarburo sería el siguiente: 
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ñ catalizador A Ñ 
Hidrocarburo > Ión carbonio + le” 


| catalizador 


Olefina + lón carbonio 


| catalizador 
Olefina + Ión carbonio 


Al final el ión carbonio puede reaccionar con otra molécula de hidrocarburo 
inicial, para dar un ión carbonio de mayor número de átomos de carbono y 
comenzar de nuevo los pasos de ruptura. Entre los catalizadores de craqueo 
empleados actualmente cabe citar las  zeolitas. Estos compuestos son 
aluminosilicatos cristalinos que contienen "cajas", a menudo de dimensiones 
moleculares, que pueden bloquear cadenas ramificadas de las moléculas y favorecer 
de esta manera la ruptura por un enlace determinado, aumentando la selectividad 


en un componente determinado. 


Los catalizadores metálicos se emplean frecuentemente en las 


deshidrogenaciones e hidrogenaciones. Un ejemplo puede ser la síntesis del NHg3_ 
que seguiría el siguiente esquema: 


Fe 
N(g) “> 2N(ad) e 
= N(ad) + H(ad)*> NH(ad) 
Hx(g) *” 2H(ad) 


NHx(2) 
| +H(ad) 


NH(ad) > NH;(ad) 
+H(ad) 
Para las reacciones de oxidación de hidrocarburos, los componentes activos 
del catalizador son óxidos metálicos (de molibdeno, vanadio, cromo, etc), capaces de 
tener varias valencias. En este caso parece que tienen lugar las siguientes 


reacciones: 


Hidrocarb. + Catalizador Oxidado => Prod. reacción + Cataliz. parcialmente 


reducido 


Cataliz. parcialmente reducido + Oy => Catalizador oxidado 


Se ha podido comprobar estas fases, marcando los átomos de oxigeno 
(isótopos) del catalizador, y comprobando posteriormente la presencia de estos 


átomos entre los productos de la reacción. 
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Se ha relacionado la actividad catalítica con el carácter conductor eléctrico 
del catalizador. Se ha comprobado que cuando el catalizador era poco conductor, no 
tenía lugar ninguna reacción, si era muy conductor la reacción más desarrollada 
era la total a óxidos de carbono y agua, y para una situación intermedia se 
formaban frecuentemente los productos oxidados deseados (ácido benzoico y 


benzaldehido a partir del tolueno). 


Para muchas reacciones, el mecanismo de reacción, que incluye la adsorción 
y desorción de los centros activos, se puede representar mediante esquemas como 


el siguiente (donde L es el centro activo): 
A+L R+L S+L 
Ty 1 Je, 
AL SRLso SL 
Este esquema podría corresponder a la reacción A«SS, donde la 


concentración de R en el fluido sería muy pequeña, por estar el equilibrio de 


adsorción muy desplazado hacia R:L 


Otras veces, con catalizadores bifuncionales, la reacción tiene lugar mediante 
la intervención de dos centros activos de diferente tipo. Por ejemplo, como se 


representa en el siguiente esquema: 

A R S 

N 1 N 

Alis RL; RL so SL 

En otras ocasiones, la molécula adsorbida puede reaccionar con otra 
molécula sin adsorber, como se muestra en el siguiente esquema: 

A e AL Cc 

| +B sl 


T— “o B:'AL S CL + D 


3.2. Teorías sobre la adsorción 


Hay dos tipos de adsorción, la adsorción física y la adsorción química o 
quimisorción, la primera se realiza por fuerzas de Van der Waals, pudiéndose 


adsorber más de una capa de moléculas sobre la superficie del sólido, mientras que 
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en la segunda la adsorción se realiza mediante enlaces covalentes químicos y 


solamente se adsorbe una capa de moléculas sobre el sólido. 


En los procesos catalíticos, previa a la formación de los compuestos 
intermedios, tiene lugar la unión de las moléculas reactivos a la superficie del 
catalizador mediante enlaces químicos. Se pueden clasificar los metales por los 


gases que son capaces de quimiadsorber. 


La clásica teoría de Langmuir, que ha justificado multitud de datos 


experimentales, se basa en las siguientes hipótesis: 
1. Los centros activos son uniformemente energéticos. 
2. Se adsorbe una sola capa. 
3. No existe interacción entre las moléculas adsorbidas. 


4. A los fenómenos de adsorción y desorción se les puede aplicar la ley de 


acción de masas. Así, para la etapa A + L «€ AL 


Si C,., es la concentración de A-L [kmol/m2], C, es la concentración en fase 
gas [kmol/m*], C, es la concentración de centros activos [kmol/m], cuando se 


alcance el equilibrio se cumplirá que: 


Ka = constante de equilibrio de quimiadsorción = a (3.1) 
A CL 


Los valores de las velocidades del proceso directo e inverso vienen dados por: 


r, (velocidad de adsorción) = k, Ca C. (3.2) 
donde k, = k¿e"Ea/RT (3.3) 
r¿ (velocidad de desorción) = ky Ca (3.4) 


donde k, = k'eEa/RT 


En el equilibrio se cumplirá: 


kC,C =k-C (3.5) 
a A L d AL 

di O (3.6) 
C,C, Ka 


donde K, es la constante de equilibrio. 
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Se puede relacionar C,, con C, de la siguiente forma: 
COFFEE? O, AFRO, O =€ M1RKC) (3.7) 
donde C: es la concentración de centros activos [centros activos/m2]. Luego queda: 


C 


=  —t 30 
E 1+K,C, (e) 
y la porción de centros activos ocupados será: 
Car Z K, Ca (3.9) 
Gs  EERyEz 


La representación de C, ,/C, en función de C, tiene la forma: 


Ca 


Figura 4. Variación de la proporción de centros activos 


C 
Si C, es alta K,+ 0,251 E =1 (3.10) 
1 
E . Car pa 
Si C, es baja K,*C, << 1 C =K,C, (3.11) 


Conviene aclarar que no siempre todos los datos experimentales se ajustan a 


esta ecuación. Se han propuesto otras ecuaciones tales como las de la: 


C 
Isoterma de Freundlich E = EC, (3.12) 


t 


Cc 
Isoterma de Temkin EA = Aln(B C,) (3.13) 


1 
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A partir del modelo de Langmuir, se pueden obtener expresiones de velocidad 
de reacción en función de la concentración de reactivos (y productos, en reacciones 


reversibles). Esto se puede ver en un ejemplo sencillo: 


En muchas reacciones catalíticas de oxidación de hidrocarburos, se ha 
comprobado que las etapas lentas, que controlan el proceso, son la quimiadsorción 
del Oz y/o la reacción química superficial entre el Oz quimiadsorbido con el 
hidrocarburo sin adsorber en fase gaseosa. Se van a deducir las ecuación cinética 


para ambos casos. El mecanismo de la reacción es el siguiente: 
1. Enlace del oxigeno a la superficie (quimiadsorción): 
O, + L > 02>L 
2. Reacción química superficial: 
OL + RH >  Compuesto-L 
3. Reacción de desorción: 
Compuesto-L —> Compuesto + L 


Si el proceso controlante es la quimiadsorción del O,, la velocidad de reacción 


vendrá dada por la etapa 1 (las demás serán muy rápidas), y será: 
r = k Co2 €. (3.14) 
Por otra parte se cumple que: 
C; =C,+CL.02 (3.15) 
Al ser la reacción química superficial más rápida Co. = O y entonces C, = C; 
Luego: 
r = kCo, C, = kC¡Co. = Kk' Cos (3.16) 


es decir la velocidad de reacción es directamente proporcional a la concentración de 


oxigeno. 


Si la etapa controlante es la reacción química superficial, la etapa 1 estará en 
equilibrio (de constante Kad), la etapa 2 será lenta e irreversible y la 3 rápida e 


irreversible. En este caso la velocidad vendrá dada por: 
r=k Co. Cr (3. 17) 
y Es =C +Co».L (3. 18) 
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Recordando que: 


C 
Ka (3.19) 
Co Cs 


y por tanto Cu E:CoCEr. 118,20) 


ya que al ser la etapa de adsorción la más lenta los reactivos y los productos 


estarán en equilibrio. Se deduce: 
C, =C,+Kg Cor Cy =C 0+ Ko Co) 5 Ch; 


Entonces 


KK Colo U ad E Es 
r=kCo., Com =KCgn Co Ci, Ko = AL RH Lo2 
+ Ko Cos IEC 


Efectivamente se han comprobado experimentalmente las dos expresiones 
deducidas. Para cualquier otro mecanismo se pueden obtener las ecuaciones 
correspondientes. Las expresiones de velocidad coincidirán con las obtenidas 
experimentalmente, siempre que sean correctas las hipótesis admitidas (mecanismo 
de Langmuir, etapa controlante, etc...), además de que las restantes etapas de la 
reacción, las fisicas, sean más rápidas. El valor de r de las ecuaciones anteriores se 
puede igualar con cualquiera de las expresiones de velocidad, basadas en la masa 
de catalizador, volumen del reactor, etc...). De ello dependerá las unidades de la 


constante o constantes. 


A partir de las ecuaciones obtenidas se puede predecir si el modelo es el 
correcto, cómo variará la velocidad de reacción con la presión total del sistema, o 


con las concentraciones de los reactivos. 


4. Ecuación de velocidad para la difusión externa como etapa limitante 


Cuando la resistencia a la difusión externa es más elevada que todas las 
demás resistencias, la difusión del reactivo desde el seno del fluido hasta la 
superficie del catalizador es la etapa limitante de la velocidad de reacción. El flujo 
de reactivo que llega hasta la superficie de catalizador viene regido por el coeficiente 
de transporte de materia k, entre el fluido y el sólido. La velocidad en este caso 


vendría dada por: 
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€. 


- [es =k (Ca, 64) [kmol / (s"m2)] (4.1) 


donde 

k¿ = coeficiente individual de transferencia de materia (((kmol / (sm2) /kmol /m3)) 
S.= Superficie externa de la partícula (m2) 

Cag= Concentración del reactivo A en el fluido (kmol/ms) 

Cas= Concentración en la superficie (kmol / m3) 


para los casos donde todas las demás etapas sean rápidas incluida la difusión 
exterior de los productos, Cas es la concentración de equilibrio del reactivo A con los 
reactivos; si las demás etapas son muy rápidas, incluida la reacción química, y si la 


reacción es irreversible, entonces Cas vale cero. 


La etapa de difusión externa de reactivo suele ser la dominante en las 
reacciones que transcurren a elevada temperatura y con partículas de catalizador 
muy pequeñas. En estos casos, la velocidad de la reacción global será la misma 
para todas las reacciones sólido-fluido, tanto para reacciones catalizadas como para 


las que no lo son. 


El problema cuando la difusión externa es la etapa controlante suele ser el 
cálculo del coeficiente de transferencia de materia, aunque se puede estimar 
mediante ecuaciones semiempiíricas donde se utilizan módulos adimensionales. Se 


puede relacionar k, con el módulo de Sherwood de la siguiente manera: 


Sh= (4.2) 


donde D= Coeficiente de difusión molecular, que se debe conocer o estimar 
mediante ecuaciones semiempíricas, d,= Diámetro efectivo de la partícula, es decir, 


el diámetro de una esfera que tenga la misma área exterior. 


Se relaciona el módulo de Sherwood con los números de Reynolds y Schmidt, 


mediante el parámetro j,; 


Sh 


4.3 
(Sc) Re pa 


Ip = 
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donde: Re = (4.4) y Sc = += (4.5) 


siendo u la velocidad superficial, calculada como caudal/sección total del reactor. 


Para lechos rellenos y para reacciones del tipo A > B en ausencia de inerte (y 


por tanto hay contradifusión) el valor de j, depende fundamentalmente del módulo 


de Reynolds, como se aprecia en la Figura 5 (véase Froment y Bischoff, 1979). 


10! 


109 
Jo es 


101 


102 


100 10! 102 103 104 
Re 


Figura 5. Transferencia de masa entre un fluido y un lecho de partículas, (relación 
jo vs Re para distintos autores, véase Froment y Bischoff, 1979). 
Existen también expresiones numéricas que relacionan j, con el módulo de 


Reynolds, como por ejemplo la ecuación de Thodos: 


A 
22" TRe P"-015) 


(4.6) 


Cuando no hay contradifusión (hay contradifusión, cuando el número de 
moles que van hacia la partícula es igual al número de moles que vienen de la 
partícula), el valor de k¿ se obtiene mediante las expresiones anteriores y dividiendo 


Por Y fa: 


Ko 


k.= (4.7) 


Y ta 


donde k¿? es el coeficiente individual calculado para contradifusión y y¡a es la 


fracción molar de A, que para un esquema de la reacción aA+bB+...> rR+sS+... 


_ (+9, y,J-(1+9, Ya.) 48) y 5, +54) lero.) 


Y fa 
dl 1+0, Ya a 
1+0,)Ya, 


(4.9) 
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En presencia de una gran cantidad de inerte y¡í se debe sustituir por un 
factor que es igual a la fracción molar de inerte que prácticamente es constante en 
la capa que rodea la partícula. En las ecuaciones anteriores, y,es la fracción molar 


de A en el fluido, y yas en la superficie. 


Con todas estas relaciones se puede estimar el valor del coeficiente k;. 


5. Tipos de difusión interna 


El caso en el que aparece una elevada resistencia a la difusión en los poros 
es un caso complicado, ya que no está relacionado de forma sencilla con las otras 


etapas. 


La expresión general de difusión de un componente A a través de una mezcla 


gaseosa A+B, viene dada por la ley de Fick: 
Ni, = Xa(N, + Ng) - CDi VXa [kmolA /m2:-s] (5.1) 


donde la densidad de flujo molar, Na, depende de dos términos, el primero de ellos 
consecuencia del movimiento general y el segundo de la difusión. Cuando no hay 
movimiento global, es decir, si existe contradifusión y ésta solo tiene lugar en la 


dirección X, se cumple que N', = -N'z y la ecuación anterior se transforma en: 


dx 
Na=-0D4o| 2) (5.2) 
dx 
y si la concentración es constante: 


dC 
N,=-=D pp ES (5.3) 


En sistemas multicomponentes, el tratamiento se complica porque el término 
de desplazamiento global no es despreciable, y las difusividades pueden depender 
de la composición. En sistemas gaseosos en presencia de una gran cantidad de 
inerte en reposo y en sistemas líquidos las ecuaciones (5.2) y (5.3) son aplicables. 
Nótese que las unidades del coeficiente de difusión, por definición, son 


M'fuido/ Mcat'S, que usualmente se indican como m2/s. 


El valor de D,¿ varía con la temperatura y la presión de diferente forma 


según el tipo de difusión. 
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Los catalizadores industriales tienen una superficie interna Sy que varía de 
10 a 2000 m2/g. Dentro de la partícula catalítica, hay una distribución de poros, 


denominándose: 
- macroporos: anchura superior a 50 nm (geles porosos, vidrios) 


- mesoporos: anchuras entre 2 y 50 nm (arcillas, tamices moleculares de 


silice-alúminas) 
- microporos: anchuras inferior a 2 nm (zeolitas, carbones activos). 


Esta distribución de poros se puede estimar mediante isotermas de adsorción 
con N2, CO», butano, porosimetría de mercurio, etc. Una distribución típica de la 
distribución de poros puede observarse en la Figura 6, donde se ha representado 
AV/Aln r en función del radio de poro "r", y siendo "v" el volumen de los poros que 
tienen un radio menor o igual a "r". En el pellet 1, la mayor parte de los poros están 
en el rango 20 a 200 Amstrong (2 a 20 nm). Se define la porosidad interna e como la 
fracción volumétrica de huecos internos. Para los catalizadores mencionados e = 
0.85 y e = 0.59. Suponiendo que los poros fueran cilíndricos, el radio medio del poro 


r se puede calcular mediante la relación: 


Z ZE 


LS 5.4 
pS, (5.4) 


siendo p, la densidad aparente del sólido (masa de sólido /volumen total incluyendo 


poros) y Sy la superficie especifica interna (m2/g). 


La difusión interna de un compuesto A a través de un catalizador puede 


tener lugar por alguno de los siguientes mecanismos, reflejados en la Figura 7. 


del tamaño de los poros, por lo que unas podrán pasar y otras no, y a la vez las que 
pasan quizás solo pueden hacerlo longitudinalmente y no transversalmente. Por 
tanto en este caso la difusividad puede variar entre cero hasta algunos valores 


finitos, pero muy bajos. Esta difusividad se le conoce con el nombre de difusividad 


incidencia como el choque molécula-molécula. Esta difusividad puede ser frecuente 


Figura 6. Distribución de tamaños de poro en catalizadores. 


D, cn?/s 


Regular 


Gases a l atm 


1 Gases a 10 atm 


Knudsen 


1 Líquidos 


Configuracional 


1 10 100 1000 10000 


rnA 


Figura 7. Tipos de difusión en procesos catalíticos gas-sólido. 


1-A través de poros muy pequeños, el tamaño de las moléculas es del orden 


configuracional. 


2-Para tamaños de poro del orden de 2 a 100 nm, el tamaño es del orden del 
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recorrido libre medio de las moléculas y el choque pared-molécula tiene tanta 
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para gases, y se conoce con el nombre de difusividad de Knudsen D,,, que se 


calcula mediante la expresión: 
ABREN 
Di (5.5) 


donde M, es el peso molecular de la sustancia que se difunde, “res el radio del poro 


donde tiene lugar la difusión y tr es el número pi. 


3- A tamaños superiores a 100 nm, tiene lugar la difusión ordinaría 


(difusión molecular por choque molécula-molécula). 


Para tamaños intermedios de poros (entre los correspondientes a la difusión 


de Knudsen y la difusión ordinaria), la difusividad D,, para el caso de AS£SB o un 


sistema con una gran proporción de inerte, se puede calcular como: 


= + (5.6) 


En los catalizadores porosos, y para un proceso A£SB, o en sistemas con gran 


proporción de inerte, la ley de difusión se aplica como 
dC 
N,=-D, (e) (5.7) 


siendo N, la densidad de flujo, definida como los kmoles de A/s que se difunden 


por unidad de área normal total, incluyendo la ocupada por los sólidos. Teniendo en 
cuenta que la sección libre de paso de los huecos es Se, el coeficiente de difusión 
sería Dyg , pero debido al aumento de recorrido que tienen que hacer las moléculas, 


se considera la siguiente relación: 


D,e 
a (5.8) 
T 


D 


donde 1 es el factor de tortuosidad (valor frecuente 3646 1/e). 


Las expresiones anteriores son válidas cuando hay contradifusión molecular 
(ASB) o el proceso tiene lugar en gran presencia de un inerte y referido 
normalmente a mezclas gaseosas. En sistemas multicomponentes debe hacerse un 
tratamiento particularizado, teniendo en cuenta que los reactivos van hacia el 
interior de las partículas, y los productos hacia fuera. El o los coeficientes de 


difusión efectivos deben de tener en cuenta estas consideraciones (consultar 
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bibliografia para ampliación de estos aspectos). Cuando hay distribución de poros, 
es posible que para los pequeños la difusividad de Knudsen sea controlante y para 
los mayores no; en este caso deben usarse las expresiones correctas (consultar 


bibliografia). 


6. Velocidad de reacción cuando la difusión interna es lenta 


6.1. Caso general de primer orden 


Se considerará el caso de una molécula de reactivo A, en un sistema de 
reacción A>B, o en un sistema que difunde a través de una pastilla porosa plana, 
en la cual se está produciendo una reacción de primer orden, como se aprecia en la 


Figura 8. 


2L 


> 
> — 


Figura 8. Pastilla porosa con entrada de reactivo A por ambas caras 


Realizando un balance de materia en régimen estacionario para el reactivo A 


en el dV =S dx de la pastilla considerado, se obtiene: 
Salida - entrada + acumulación = generación — [kmolA/6s] 


Nal... S-N,1,S+ 0 = r,:dV = r,:S:dx (6.1) 


x+dx 


siendo Na la densidad de flujo de A (kmol A/(s.m2)) y ra la velocidad de reacción de A 
(kmol A/(s m3 cat)), que se va a suponer que es de primer orden, por lo que se 


deduce: 


d| —D, ee 
IN, - 
po =r, =-k-C, (6.2) y de aquí a (6.3) 
donde k es la constante de la reacción de primer orden [m3fuido/ (m%cat.s)], D, es la 
difusividad a través de la mezcla gaseosa (como ya se ha indicado, con unidades de 


M3fñuido / Meat"S) . 
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Por lo tanto, en este caso particular: 


d?*C A 
=D, a =-kC, (6.4) o  D, — =kC,=0 (6.5) 
da Tk Ñ 
A C4 20 (6.6) 


La solución general de esta ecuación diferencial es: 
Ca= M;¡e mx + Moe nx (6.7) 


donde el valor de m viene dado por: 


e 
m[ ] (6.8) 


Este módulo m, según lo expuesto anteriormente, tiene las unidades de 
1 /Moat. 
Los valores M, y M) son unas constantes que pueden ser determinadas por 


las condiciones de contorno. Una de estas condiciones es que en la superficie de la 
pastilla la concentración de A es igual a la que existe en el seno del fluido. La otra 
supone que la variación de la concentración en la mitad de la pastilla porosa es 


cero. 
1) C,=Ca) parax=0 (6.9) 
2) dC,/dx = O, para x = L (6.10) 


Si tenemos en cuenta la simetría de la pastilla, la condición 2 no es más que 
indicar que en el centro la concentración de A presentará un minimo. La condición 


1, es también lógica. 
A partir de aquí, se obtiene: 


Ex pont Cs ent 
M3 == Ms ===. (611) 
+ 


y por lo tanto: 


C, RS SL a Ñ cosh[m(L — x)] 


6.12 
Cá enter cosh(mL) e 
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(Función matemática: cosh a=(e%+e-0)/2, senh a=(e%-e-0)/2; tanh a=(ed-e-0) / (¿4-01 


A la vista de esta ecuación se puede deducir que se produce un descenso de 
la concentración de reactivo conforme éste va avanzando hacia el interior de la 
pastilla. Este descenso depende de una magnitud adimensional mL, que se define 


como módulo P de Thiele. En la Figura 9 se representan distintas gráficas C,/C,, vs 


(1-x/L) para diferentes valores del módulo de Thiele. 


Si solo existe difusión a través del eje x, se pueden calcular las densidades de 


flujo de entrada y salida y, multiplicando por la sección, los correspondientes flujos. 


Puede observarse que cuanto mayor es el módulo de Thiele, o sea, cuanto 
más rápida sea la reacción frente a la difusión, más pronunciada será la 
disminución de la concentración de A. Sin embargo, para valores del módulo de 
Thiele pequeños, la difusión a través de la pastilla es muy rápida y la concentración 


de A no varía mucho. 


Ñ 0.8 0.6 0.4 0.2 Í 0 
Centro de la partícula 


Superficie exterior 1-(x/L) 


Figura 9. Variación de la concentración de A frente a la distancia x. 


Para determinar la disminución de la velocidad de la reacción en el caso en el 
que la difusión a través de la pastilla porosa sea importante, se introduce una 


nueva magnitud, el factor de eficacia, que se expresa por: 


generación real (kmol A/ s) velocidad real media (kmol A((sm” cat)) 


generación sin efecto de difusión interna(kmol Als) velocidad media sin efecto de difusióninterna(kmol A((s m* cat)) 
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v 


[epa , 
E A J (=riddL J (—D, aca ML +D, 23 
Ya V 0 0 dx X |,=0 


= A (0318) 


¿ds Er; ñas A 07 KC ,.L KC, L 


donde ls r;) y E r,) se expresan en (kmoles de A/(m'.as), V es el volumen de la 


partícula catalítica y Sita representa la superficie total por la que hay entrada de 


reactivo A (en el caso considerado Sita = 2 S). Nótese la diferencia entre la velocidad 
real media (- r,) que puede ser determinada experimentalmente, y la velocidad 


Ll r¡Jen cada capa de la partícula porosa. 


Luego Er,)=kC,, 7 (6.14) 


lo cual es una expresión muy cómoda para trabajar, ya que el factor de eficacia 


modifica la velocidad de reacción según el efecto de la difusión. 


Tomando la ecuación del perfil de concentración obtenido (ec. 6.12) se obtiene: 


dC, _ Cas 
dx cosh(mL) mpserhln a) 


que para x=0 proporciona: 


dC, _  mC,, 


senh(mLÑ) ==mC ,. tanh(mL 
dx |, cosh(mL) ( ) E Sid 


Introduciendo esta ecuación en (6.14) obtenemos: 


tanh (mL) 


E (6.15) 


Teniendo en cuenta los datos de esta gráfica (6.2) se pueden hacer dos 


aproximaciones para simplificar la velocidad de reacción obtenida: 


- Cuando no hay resistencia a la difusión en la pastilla, (P=mL<0.4), nl, y: 
=r,)=*C,, (6.16) 
- Cuando la resistencia a la difusión es importante, (Pb=mL>4), n=1/mL, y: 


Er Es la D, ya] 
mL L 


Er)- 


(6.17) 
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En principio, podría parecer que la difusión en la pastilla puede considerarse 
como un factor multiplicativo correctivo de cualquiera de los casos anteriores, 
utilizando n. Pero esto no es así, ya que este factor lleva implícito en su coeficiente 
cinético un término de reacción en la superficie además del término de difusión en 


la misma. 


El problema planteado para la pastilla porosa con transporte de A a través de 
las dos caras, se puede plantear para otros casos: pastilla porosa permeable a los 
reactivos y productos por una sola cara, pellets cilindricos y partículas esféricas. 
Para el caso de la pastilla porosa permeable sola por una cara, la solución es 
exactamente la misma que el problema resuelto anteriormente pero considerando 
"L" al espesor total. Para los otros dos casos, pellets cilindricos y partículas 
esféricas, existen soluciones analíticas exactas. No obstante, introduciendo una 
longitud característica L (definida como L = Volumen de la partícula/Superficie 
externa) las soluciones analíticas son muy próximas, tal como se indica en la Figura 
10 (en esta figura, L es la longitud característica, no es el radio). Esta longitud 


característica se define por tanto como: 


- Placa plana sin efectos de extremos con difusión por ambas caras: 
L = (mitad del espesor) 
- Placa plana sin efectos de extremos con difusión por una cara: 
L = (espesor) 
- Pastillas cilíndricas largas: L = radio /2 
- Pastillas esféricas: L = radio/3 
- Para pastillas que no se puedan clasificar dentro de las anteriores: 
L = (volumen de la pastilla) / (superficie exterior disponible para la entrada y 


difusión del reactivo). 


1.0 —===="=x == 
n 0.8 


06 ——— Esfera 


O Cilindro 
ESDEZES Lámina 
0.4 


0.2 


0.1 


0.1 1 10 


Figura 10. Variación del factor de eficacia frente al módulo de Thiele 
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El error más grande se obtiene en las zonas centrales de la curva, donde se 


comprobó que está por debajo del 18%. 


6.2. Caso en que el producto tenga problemas de difusión interna 


Otro caso posible, poco analizado en los libros de texto, es el del producto B 
(siguiendo con el esquema A>B) con una resistencia no despreciable a la difusión 
dentro de la pastilla porosa. Se analiza a continuación el caso de una pastilla 
porosa plana, con salida de producto B por ambas caras y en la que se está 
produciendo una reacción de primer orden. Realizando un balance de materia en 
régimen estacionario para el producto A en el dV =S dx de la pastilla considerado, 


se obtiene: 
Salida - entrada + acumulación = generación — [kmolB/s] 


Noli. S+N¿L,S+ 0 = rj:dV = r,"S:dx (6.18) 


x+dx 


En este caso el signo de los flujos molares de B es el contrario que el de los 
Na. Si, como anteriormente, se trata de una ecuación de primer orden, y utilizando 


la ley de Fick para la difusión de B: 


ry=KC as ya que Ca es constante dentro de la pastilla 


ley de Fick para el componente B 


(6.19) 


Como vemos, en este caso la derivada segunda de la concentración es una 


constante. Si se integra la ecuación, resulta: 


dC, kC 
Ma Nas (6.20) 
dx D, 
KC, x 
B ES 0, (6.2 1) 


e 


donde 61 y 62 son constantes de integración que es necesario determinar. 
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Las condiciones de contorno del sistema son similares a las condiciones 
estudiadas en el caso del reactivo A con problemas de difusión. Asi, se tiene: 

1) Cá= Cay parax=0 (6.22) 

2) dCs/dx = 0, para x = L (6.23) 


En este caso la condición 2 indica que en el centro de la pastilla, la 


concentración de B presentará un máximo. La condición 1, es también lógica. 


Aplicando las condiciones a la ecuación (6.21) se llega a: 
==> y 0,=C>, (6.24) 


con lo que, sustituyendo y reordenando: 


KC, (í 6 : 
us E e al 008 
Da Ca  DL2 


Ss 


e 


Esta ecuación se ha representado en la Figura 11 para varios valores de 


m=(k/D.)!/2, suponiendo Cgs=0.1 y Cas=1. 


1 Ñ 0.8 0.6 0.4 0.2 Po 


Superficie exterior 1-(x/L) Centro de la partícula 


Figura 11. Variación de la concentración de B frente a la distancia. 


Como era de esperar, cuanto mayor es el módulo de Thiele (mayor velocidad 
de reacción química frente a la difusión), los perfiles de concentración se hacen más 


pronunciados. 
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6.3. Mezclas de partículas de formas y tamaños diferentes 


Cuando se tiene una gran cantidad de partículas de formas y tamaños 
distintos, se debe aplicar el factor de eficacia medio que se puede calcular de la 


siguiente manera: 
rd a AN (6.26) 


donde f,, f,,... son las fracciones volumétricas de partículas con tamaños distintos 


respectivamente. 


6.4. Cinética de reacción arbitraria 


Para una reacción reversible de primer orden, con difusión rápida de los 
productos y por tanto con concentración constante dentro de la partícula, se puede 


deducir que el factor de eficacia viene dado por 


Y 


mL = ll (6.27) 
D, Cae 


siendo £x. el grado de conversión de A en equilibrio con los productos de reacción 
(en el caso de que la difusión de los productos sea también importante, se deben 


plantear las correspondientes ecuaciones y resolverlas numéricamente). 


Aris y Thiele construyeron una gráfica que representaba el factor de eficacia, 
n, frente al módulo de Thiele, mL, para cinéticas de orden distinto a uno. Se 
encontró que para reacciones con orden n, estas curvas se aproximaban mucho a la 
curva correspondiente de primer orden, (Figura 12), generalizando el módulo de 
Thiele mediante la expresión: 

n-1 % 
mL=L al (6.28) 

2D, 


siendo Cas la concentración de Ca en la superficie exterior de la partícula. 


En ambas expresiones (6.27) y (6.28), D,¿ es la difusividad efectiva del fluido 


en el sólido. 


DISEÑO DE REACTORES HETEROGÉNOS 27 


Figura 12. Factor de eficacia y módulo de Thiele generalizado. 


6.5. Ensayo de un catalizador con efectos importantes de difusión interna 


Supóngase una reacción en la que se ha determinado la velocidad 
experimental y en la que la resistencia a la difusión externa es despreciable. En este 
apartado se intentará determinar si la difusión interna puede o no ser importante. 


Para ello, se estudian por separado las reacciones de orden uno y las de orden n. 


1. Reacciones de primer orden.- La velocidad de reacción cuando se 


considera la difusión en una partícula catalítica es: 
=r,)=kC,, y (6.29) 


donde las unidades de (-r,) dependen de las unidades en que se exprese la 


constante cinética. 


Teniendo en cuenta el módulo de Thiele, se puede deducir que 
47 =p (mLY (6.30) 


Este es el llamado módulo de Weisz (We). 


A 0 (6.31) 
Do. 
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gb 0 


2. Reacciones de orden n.- Cuando la ecuación cinética es del tipo 
ra=k-f(Cas), se puede introducir la expresión generalizada del módulo de Thiele 


descrita anteriormente y se obtienen las expresiones correspondientes, del tipo: 


2 
=Y, E 


D, Ca 


=7] (mL) (6.33) 


7. Efecto de la temperatura en la velocidad de reacción 


Cuando la reacción es tan rápida que no permite que el calor absorbido o 
desprendido en las proximidades de la superficie del catalizador o dentro de éste se 
disipe, habrá que tener en cuenta la contribución de este efecto en la velocidad de 


reacción. 


Si la reacción es exotérmica convendrá favorecer las condiciones para que se 
produzca una reacción no isotérmica, ya que se tendrán temperaturas elevadas 
cerca de la superficie del catalizador y la reacción se llevará a cabo más 
rápidamente, (siempre que el aumento de la temperatura no conlleve a condiciones 
adversas como la posible descomposición del reactivo). Para reacciones 
endotérmicas, siguiendo el mismo razonamiento, convendrá trabajar siempre en 


condiciones isotérmicas. 
Se pueden considerar dos gradientes de temperatura distintos: 
7. 1. Gradiente externo de temperatura entre el medio y la superficie de la 
partícula 


Se presenta cuando existe una diferencia de temperatura entre la superficie 
del catalizador y el fluido externo. El calor cedido al fluido debe ser igual al 


generado en la superficie del catalizador: 


Oseneraso Va (Eta [E AB,) (7.1) 


Obecióo = NS a (Y, —T,) (7.2) 
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donde Q viene dado en J/s, (-r"») en kmol A/(s m? catal), (AH) la entalpia de 
reacción en J/kmol A, Voat el volumen de la pastilla, Scar la superficie exterior de la 
pastilla, Ts la temperatura de la superficie de la partícula, Ty la temperatura del 
fluido que rodea a la partícula y h el coeficiente individual de transporte de calor 
(J/(sm2*C)). Igualando estas dos expresiones, se puede llegar a: 


AT 


película 


y 2 Er JEa4,) 
h 


=(r, —T, (7.3) 


siendo L la longitud característica Voeat/ Scar. 


El valor de h puede calcularse en algunos casos introduciendo el parámetro 


ju: 


j -( L pre (7.4) 
HT . 
C,up 


siendo (u:p) el gasto másico, Pr el número de Prandt = EA , C, el calor específico 
K 


term 


del fluido, y Kterm la conductividad térmica del fluido. 


En la Figura 13 se muestran los resultados obtenidos por distintos autores al 


correlacionar j,, frente al número de Reynolds para un lecho relleno. 


100 


E y 
Ju 10 === == 
ap 


SS 


102 


10 10! 102 103 104 
Re 
Figura 13. Representación de ju frente al número de Reynolds. Curva 1: Wilke 
y Hougen. Curva 2: Baumeister y Bennett (a) para dt/dp>20, (b) para dt/dp<20. 
Curva 3: Glaser y Thodos. Curva 4: Acetis y Thodos. Curva 5: Sen Gupta y Thodos. 
Curva 6: Handley y Heggs. 
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Con los datos de la gráfica anterior puede ser determinado j,, y a partir de él, 


el valor del coeficiente individual, h. Para reacciones de primer orden en una 


reacción entre gases, en el caso en que exista contradifusión se cumple que: 


armo e Jo, Mus) 
M,C» 

En esta ecuación Mm es el peso molecular medio del gas. En el caso en que 
no exista una diferencia notable entre las fracciones molares de A en la superficie y 
en el seno del fluido, esto es, en el caso en que no existan problemas de difusión 
externa, la diferencia de temperatura también sería cero. Si no existe contradifusión 


la ecuación correspondiente vendría dividida por el factor yya. 


La máxima diferencia de temperatura se producirá cuando la reacción sea 
completa y rápida (esto es, yas=0), por lo que: 
E ao! 


(AT)... = 07 
M,C 


mp 


que en el caso en que no exista contradifusión se modificaría a: 


(ar) AAA 


ds Ya m C, M, C, Ora 


donde %, es la expansión molar del compuesto, es decir, la diferencia de la suma de 


los coeficientes estequiométricos de reactivos y productos, y T la presión total. 


7.2. Gradiente de temperatura en el interior de la partícula 


Se produce cuando existen diferentes temperaturas en distintos puntos de la 
pastilla de catalizador. Prater demostró que este efecto es el mismo para distintas 
formas geométricas de partícula y para reacciones de diferentes órdenes. El estudio 
de este fenómeno se basa en que la temperatura y la concentración dentro de la 
partícula vienen determinadas por ecuaciones diferenciales semejantes. Prater 


demostró que para cualquier punto de la pastilla se cumple: 


Alano can, 05 


dx dx 


siendo k,, la conductividad efectiva del catalizador, que integrando: 
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AT Tr ) E [D, (C,, == nta JE AH, )] 


partícula al e A 
k 


(7.6) 
ef 


De nuevo se obtiene que en el caso de no resistencia a la difusión, es decir, 
cuando Cas=Cacentro, la diferencia de temperatura es cero. Esta expresión permite 
conocer el valor de ATparticula si se conoce Teentro y Un valor limite cuando Cacentro=0. 
Para el gradiente de temperatura dentro de la partícula, se ha representado el valor 
del factor de eficacia n, frente al módulo adimensional de Thiele generalizado, para 
reacciones de primer orden obtenido por simulación numérica (Figura 14). En esta 
figura, la curva correspondiente a pP=0 corresponderá con la representada 
anteriormente en el caso del régimen isotermo (figuras 10 y 12). Puede observarse 
que para reacciones exotérmicas, el factor de eficacia global puede ser mayor que la 
unidad. Este comportamiento está basado en que con el suficiente incremento de 
temperatura causado por las limitaciones en la trasferencia de calor, el aumento en 
la constante cinética compensa la disminución en la concentración de reactivos, de 


forma que la velocidad de reacción es mayor que en las condiciones de la superficie. 


Se puede observar que, para un intervalo estrecho de valores del módulo de 
Thiele, se produce así la multiplicidad de estados estacionarios. Se puede deducir 
que los dos estados extremos son estables y es inestable el estado intermedio. En la 
Figura 14 se ha presentado los casos A, B y C para un determinado valor del 
módulo de Thiele, es decir, para una reacción dada, con un catalizador dado, de los 
cuales sólo A y C son estables. Cuál de los dos prevalezca es función de por dónde 
se llegue al sistema y de las condiciones iniciales. En la práctica ordinaria es 
extraño encontrar estos casos. También es importante indicar que en ciertas 
regiones cercanas a la parte vertical de la curva, un pequeño cambio en el módulo 


de Thiele puede causar un gran cambio en el factor de eficacia. 
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0.1 


0.01 


0.001 1 


Figura 14. Factor de eficacia en procesos no isotermos para sistema con ley 
cinética de primer orden. B=ATmax/ Ts, y=E/RTs. Eje X: Modulo de Thiele en la 


superficie. Eje Y: Factor de eficacia. 


En la bibliografia, se pueden localizar muchas otras expresiones exactas oO 
aproximadas para calcular el factor de eficacia en diferentes casos: sistemas 
distintos de ASR, cinéticas de orden diferente a la unidad, control de difusión 


interna y externa de materia simultáneamente, etc... 


8. Tipos de reactores 


El contacto de un reactivo gaseoso con un catalizador sólido se puede hacer 
de diversas formas, cada una de las cuales posee ventajas y desventajas sobre el 


resto. 


Según el tipo de contacto los reactores pueden dividirse en dos grandes 


grupos: 


1. Reactores de lecho fijo. 
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2. Reactores de lecho fluidizado. 
Si se comparan entre sí, se puede observar lo siguiente: 


> La circulación de gases a través de lechos fijos se aproxima al flujo en 
pistón, y es completamente diferente a la circulación en lechos fluidizados, donde el 
flujo no está bien definido pero sí se sabe que es muy diferente al flujo pistón y que 
puede presentar un gran bypass, debido a que parte del gas pasa en forma de 
burbujas. Esto hace que el contacto no sea efectivo necesitándose mucho 
catalizador si se quieren conversiones altas. Si para el caso estudiado la eficacia del 


contacto es muy importante es mejor usar un reactor de lecho fijo. 


> Es difícil controlar la temperatura en los lechos fijos sobre todo si éstos 
son grandes ya que su conductividad de calor es baja, por lo tanto, en reacciones 
muy exotérmicas se formarán zonas calientes o frentes calientes móviles que 
pueden perjudicar al catalizador e incluso desactivarlo. Esto no ocurre en los lechos 
fluidizados donde al existir una mezcla rápida se puede igualar fácilmente la 
temperatura trabajando prácticamente en condiciones isotermas. Si el problema en 
cuestión requiere un alto control de la temperatura no pudiendo ésta variar mucho, 


lo mejor es usar un reactor de lecho fluidizado. 


> En los lechos fijos no se puede emplear un catalizador de un tamaño muy 
pequeño ya que se formarían tapones y habría grandes pérdidas de presión; sin 
embargo en lechos fluidizados esto no ocurre y se pueden emplear catalizadores de 
pequeño tamaño, lo que hace que se utilicen los mismos de un forma mucho más 
eficaz. Por lo tanto si la reacción química es muy rápida y las etapas lentas son la 
difusión externa e interna, interesará que el catalizador sea lo más pequeño posible 
para que exista la mayor superficie de reacción disponible y por tanto se debe usar 


lecho fluidizado. 


> Si el catalizador utilizado en el proceso se desactiva fácilmente y hay que 
regenerarlo muy a menudo, es mejor usar lecho fluidizado ya que al comportarse 
los sólidos como un fluido se pueden manejar mucho más fácilmente y puede haber 


una recirculación de sólidos desde un lecho de reacción a un lecho de regeneración. 
Los distintos tipos de diseño de estos dos grupos de reactores principales son: 
a) Para lecho fijo 


1. Lecho de relleno típico con refrigeración externa, lo que hace que el 
interior del reactor esté mucho más caliente que las paredes. Si en su interior se 
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produce una reacción exotérmica, se producen gradientes radiales de todos los 
tipos (de materia.- al producirse la reacción más rápidamente en el interior por 


estar más caliente, de calor.- por la diferencia de temperatura). 


2. Lecho de relleno pero donde se ha introducido también refrigeración 
en el interior del lecho para evitar el problema de las zonas calientes que se ha 


visto en el punto anterior. 


3. Lecho de relleno dividido en distintas etapas y donde se produce 
refrigeración entre cada una de ellas, pudiendo estar todo esto colocado en la 
misma carcasa o instalado por separado. Como la velocidad es mucho más rápida 
al principio se suele colocar una menor cantidad de catalizador en la primera etapa 
e ir aumentando esta cantidad progresivamente en las distintas etapas, para evitar 


una alta conversión en la primera etapa y un aumento muy grande de temperatura. 
ii) Para lecho fluidizado: 


1. Lecho fluidizado simple usado cuando el catalizador es estable y no 
necesita regenerarse. Existe un intercambiador de calor que se introduce en el 
reactor para un mejor control de la temperatura y evitar un gran calentamiento o 


enfriamiento del mismo. 


2. Lecho fluidizado usado cuando el catalizador debe regenerarse 
continuamente. En este caso se suelen utilizar dos lechos fluidizados: uno de 


reacción donde se introducen los reactivos, y otro de regeneración del catalizador. 


3. Lecho fluidizado de múltiples etapas con circulación de sólidos y 
mezcla reactiva en contracorriente, donde se evitan aumentos grandes de 


temperatura y se va a regenerar el catalizador. 


4. Reactor de lecho móvil. Un reactor intermedio entre los dos grandes 
grupos anteriores (lecho fijo y lecho fluidizado) es el reactor de lecho móvil, parecido 
a los de lecho fijo en que el flujo del sólidos es como en el flujo pistón. No se puede 
usar catalizadores de tamaño pequeño (para evitar apelmazamientos). Por otra parte 


permite un fácil manejo del catalizador con lo que el coste se reduce. 


Para un buen diseño hay que tener en cuenta numerosos factores, pudiendo 
ocurrir que el mejor diseño se obtiene al combinar alguno de los tipos anteriores. 
Ejemplo: si se quisiera conseguir una alta conversión en reacciones muy 


exotérmicas lo mejor será usar primero un reactor de lecho fluidizado con lo que se 
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consigue una homogenización de la temperatura y después ya un lecho fijo con lo 


que se aumenta la conversión ( puntos 1 y 2 del comienzo de la introducción ). 


Las dificultades del diseño de los reactores catalíticos recaen en como tener 
en cuenta el comportamiento no isotérmico de los lechos fijos y el flujo no ideal de 


los lechos fluidizados. 


9. Diseño de reactores de lecho fijo con relleno 


A continuación se va a abordar el diseño de un lecho fijo, que es una de las 
configuraciones más simples y frecuentes. En primera aproximación, se supone que 
el flujo de fluido a través de un lecho de relleno es próximo al de flujo pistón. El 


balance de energía se puede escribir como: 
Salida - Entrada + Acumulación = Generación 
Y ln, », +d(n, n,) [En h, +40 |+o=0 (9.1) 
j=l j=1 


o bien: 


d En ,,)-40=0 (9.2) 


donde ny es el flujo molar y hj es la entalpia molar del componente j. En función de 


la conversión molar (X;, kmol/s) tenemos: 
R 
n,=n, + 0,X', (9.3) 
i=l 
Sustituyendo la ec. (9.3) en la (9.2): 


dE (1.*Ea,x,)a,|-d0-=0 (9.4) 


j=l ¡=l 


Ss R s 
2 + d Ex, Eo, »,)-49=0 (9.5) 
JE iS J> 


Considerando despreciable la entalpia de mezcla, se acepta que: 


dh ¡ 


=p, (9.6) 


donde Cp; es el valor molar. Además la entalpía de reacción es: 
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Y 0h, =AH, (9.7) 


j=1 


Al sustituir queda: 


Ss R 
Ca +d Ex, a, )-40=0 (9.8) 


j=1 ¡=l 
s R R 

dT Y'n,C, +9 X' d(AH,)+ AH dX''-dQ=0 (9.9) 
j=1 ¡=1 ¡=1 


Si se tiene en cuenta que: 


d(AH,)=AC, dT= (Lac, la (9.10) 


T 
AH, = AH] + | AC,, dT (9.11) 
la ecuación queda como: 


PR YX, AC, ara +fac, ar ja -dQ=0 (9.12) 


j=l i=l i=l 


Si solo existe una reacción, el sistema es adiabático, y ACp¡=0, queda que: 


aT Y n, C, +AH,dX'=0 (9.13) 
j=l 
de donde: 
pa (9.14) 


donde AT es la diferencia de temperatura entre la corriente de salida y la de entrada 
[KK], e AH, es la entalpía de reacción por kmol de reactivo j [kJ/kmol] a la 
temperatura de la reacción, n,, es el caudal molar inicial del reactivo j [kmol j/s]. En 


función del grado de conversión del reactivo limitante sería: 


AT As El AH, MiS ' 


=- o (9.15) 
HO A 
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Lo que predice una recta (denominada línea de operación adiabática) en el 


espacio AT vs. Ex. 


Supongamos ahora que el reactor es de tipo tanque agitado. Esto sería un 
recipiente relleno de catalizador sólido que se está alimentando continuamente con 
un fluido que reacciona, a la vez que el sistema está en agitación constante. El 


balance de energía en este caso seria: 


Salida - Entrada + Acumulación = Generación 
Y n,h,- Y nyh,)+d0+0=0 (9.16) 
j=l j=l 


que para un sistema adiabático con una sola reacción en la que AC;;¡=0 se llega a: 


AH, 


AT = x' (9.17) 


ecuación que es idéntica a la obtenida en el caso de reactor de flujo pistón. La única 
diferencia es que en el reactor de tanque se llega directamente al estado final que 
predice la ecuación, mientras que en el reactor de flujo pistón se pasa por todos y 


cada uno de los puntos que van desde el estado inicial al final. 


Estos reactores adiabáticos de lecho de relleno por etapas son sistemas 
versátiles que se pueden aproximar a cualquier progresión de temperatura deseada. 
El cálculo y diseño de estos sistemas es sencillo y se puede esperar que el 


funcionamiento real se aproxime mucho al predicho. 


9.1 Temperatura óptima 


Se podría considerar como óptima la temperatura que hace mínimo el valor 


de V/n,, para una determinada conversión del reactivo en un reactor en continuo. 


Para cada sistema reaccionante siempre hay una temperatura a la cual la velocidad 
de reacción es máxima. Antes de continuar es conveniente aclarar los efectos que 
tienen la presión y la temperatura en las reacciones. A partir de la termodinámica 


se puede deducir las siguientes conclusiones (ilustradas en parte en la Figura 15): 


1. La constante de equilibrio termodinámica no depende de la presión del sistema, 
de la presencia o ausencia de inertes, ni de la cinética de la reacción, pero depende 


de la temperatura del sistema. 
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2. Aunque la constante de equilibrio termodinámica no depende de la presión o 
presencia de inertes, la concentración de equilibrio de los compuestos y la 


conversión de equilibrio de los reactivos pueden verse influidas por estas variables. 


3. K>>1 indica que la conversión puede ser prácticamente completa y la reacción 
puede considerarse irreversible. K <<1 indica que la reacción no transcurrirá en 


una extensión apreciable. 


4. Para un aumento de temperatura, la conversión de equilibrio aumenta para las 


reacciones endotérmicas y disminuye para las exotérmicas. 


5. Para un aumento de presión en las reacciones en fase gaseosa, la conversión 
aumenta cuando el número de moles disminuye con la reacción, y la conversión 


disminuye si el número de moles aumenta con la reacción. 


6. La disminución de inertes para todas las reacciones actúa del mismo modo que 


un aumento de presión en las reacciones en fase gaseosa. 
Era 
LO A 


In K>>1 — casi irreversible 
/ tt 


> K<<l E sólo es posible 
conversión muy baja 


Figura 15. Efecto de la temperatura en las reacciones. 


En reacciones irreversibles la velocidad aumenta siempre con la temperatura 
para cualquier composición, por consiguiente la velocidad máxima corresponde a la 
temperatura máxima permisible (limitada por los materiales de construcción, 
aparición de reacciones secundarias, etc...). En reacciones reversibles endotérmicas 
una elevación de la temperatura hace aumentar tanto la conversión de equilibrio 
como la velocidad de reacción, y por tanto en este caso, también hay que emplear la 
temperatura más alta permisible. En las reacciones reversibles exotérmicas 
aparecen dos factores opuestos: al elevar la temperatura aumenta la velocidad de la 


reacción directa pero disminuye la conversión máxima de acuerdo con la variación 
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de la constante de equilibrio; la progresión óptima corresponde a una variación de 
la temperatura empezando por una temperatura alta que va disminuyendo al 


aumentar la conversión. 


Considerando la expresión del balance de energía para un sistema adiabático 
(ec 9.14), se tiene que la variación de la conversión o del grado de avance frente a la 


temperatura debe ser prácticamente lineal, tal como se muestra en la Figura 16. 


Condiciones isotermas 


> 


Conversión completa 
del reactivo limitante 


Más 
inertes 
>; 


4 Más inertes 
E t—, —>/ Menos inertes 
: Menos E 

4 inertes 


Tent= temperatura alimento 


T 
Figura 16. Variación del grado de conversión vs. temperatura. 


En la figura se muestra la forma de las curvas del balance de energía, tanto 
para las reacciones endotérmicas como para las exotérmicas y tanto para reactores 
de mezcla completa como de flujo pistón. Estas curvas pueden ser casi rectas en 
muchas ocasiones. Para el reactor de flujo pistón la temperatura del fluido en el 
reactor se desplaza a lo largo de la curva; para flujo en mezcla completa, el fluido 
alcanza inmediatamente el valor final que se lee sobre la curva. Estas lineas se 
denominan líneas de operación adiabática del reactor. En estas curvas al 


aumentar los inertes aumenta C, y se acercan más a la perpendicular, al eje de 


abscisas. Hay que conocer, que una línea totalmente perpendicular al eje indica que 


la temperatura no varía durante la reacción (reacción a temperatura constante). 


El tamaño del reactor necesario para un fin determinado se calcula del modo 
siguiente: para flujo en pistón se tabulan los valores de la velocidad para distintos 
valores de €, leidos sobre la linea de operación adiabática, se representa 1/(-ra) 


frente a ¿a y se integra; para flujo en mezcla completa se utiliza simplemente la 


O Juan A. Conesa. Universidad de Alicante 


REACTORES CATALITICOS 40 


velocidad correspondiente a las condiciones existentes en el reactor. Este 
procedimiento resuelve a la vez los balances de materia y energía. Desplazándonos 
a lo largo de la linea de operación adiabática se satisface el balance de energía; 
calculando las velocidades a lo largo de esta linea y evaluando V/n». se satisface el 


de materia. 


El mejor método de operación en condiciones adiabáticas para un reactor de 
flujo pistón se encuentra trasladando la línea de operación adiabática (variando la 
temperatura de entrada), de modo que la velocidad alcance el valor medio más 
elevado, tal como se observa en la Figura 17. En esta figura, el lugar geométrico de 
velocidades máximas se obtiene representando curvas de velocidad constante y 
uniendo los puntos máximos. Para calcular las curvas de velocidad constante se 
dan valores a la velocidad de reacción y temperatura, y con la ley de velocidad se 
calcula la conversión, buscando el lugar geométrico de los máximos. En este lugar 


geométrico la velocidad es mayor cuanto mayor es la temperatura. 


FLUJO PISTÓN . MEZCLA COMPLETA 


—— Línea de balance 


Atos, Lugar geométrico de 
velocidades máximas 


as Conversión deseada 


Temperatura de entrada 
T¡'=demasiado baja 
T,”=correcta 
T,””=demasiado alta 


Ep TI, MEN 


Figura 17. Localización de la línea óptima de operación adiabática. En flujo 


pistón se necesitan tanteos previos, mientras que en mezcla completa no. 


Para las reacciones endotérmicas esto implica que se debe partir del valor 
más alto permisible de la temperatura; para reacciones exotérmicas, se debe operar 
a ambos lados del lugar geométrico de velocidades máximas, y mediante unos 
cuantos ensayos se determina la temperatura óptima de entrada que hace minimo 


el valor de V/n,,. Para un reactor de mezcla completa se ha de operar en el lugar 


geométrico de las velocidades máximas. 
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Cuando la velocidad disminuye progresivamente con la conversión se ha de 
utilizar flujo en pistón; éste es el caso de las reacciones endotérmicas. Para las 
reacciones exotérmicas que experimentan una gran elevación de la temperatura, la 
velocidad aumenta desde un valor bajo hasta un valor máximo que corresponde a 


una conversión intermedia £,, y después desciende; para este último caso, son 


mejores los dispositivos con recirculación o alimentación en varias etapas. 


A conclusiones similares se puede llegar usando la ecuación del balance de 
energía. Para valores pequeños de C,/-AH, es mejor usar mezcla completa. Para 


valores elevados de C¿/-AH, es mejor flujo pistón, tal como se observa en la Figura 
18. 


Recta de operación, 
pendiente = -C,/CAH,) 


Pequeña pendiente, conviene 
utilizar mezcla completa 


Pendiente grande = la velocidad disminuye 
continuamente, es mejor flujo pistón 


Figura 18. Tipos de reactores que conviene usar en cada caso 


Como los sistemas no suelen ser adiabáticos, para que se aproximen lo 
máximo a las condiciones ideales será necesario suministrar o absorber 
deliberadamente calor del reactor y también hay que considerar las pérdidas de 


calor a los alrededores. 


Las condiciones adiabáticas se pueden conseguir trabajando con múltiples 
etapas, cada una de las cuales se considera adiabática, con enfriamiento entre las 


mismas. 


Para cualquier número de etapas, las condiciones óptimas de operación se 


reducen a hacer mínima la cantidad total necesaria de catalizador. Se pueden 
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presentar diferentes casos, donde la velocidad de reacción en los lechos de rellenos 


se desarrollan entorno a la zona de velocidad máxima. En la Figura 19 se presentan 
diferentes casos. 


Figura 19. Diferentes alternativas para reducir la masa de catalizador y/o 


conseguir el grado de conversión deseado. 


